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L'invention est relative a un precede de lissage d*une propridte de subsurface 
dans une structure geologique representee par des mesures sismiques. 

L'invention est egaiement relative h un dispositif de lissage d'une propriety de 
subsurface dans une structure geologique representee par des mesures sismiques. 

5 Uinvention est enfin relative k un produit-programme d ! ordinateur permettant 

le fonctionnement d'un dispositif programmable de lissage d'une propriety de 
subsurface dans une structure geologique representee par des mesures sismiques. 

Le document WO 01/63323 Al d£crit un precede de traitement de donnees 
sismiques comprenant les etapes consistant h : obtenir un volume de donnees 

10 sismiques couvrant un volume predetermine de terre ; determiner pour chaque voxel 
du volume de donnees sismiques Torientation locale des donnees sismiques ; 
determiner pour chaque voxel s'il existe un bord dans le voisinage, et exdcuter une 
operation de lissage sur chaque voxel dans le volume de donnees sismiques, dans 
lequel la direction de Toperation de lissage est Torientation locale des donnees, et 

15 dans lequel l'operation de lissage ne depasse pas le bord, de mantere k obtenir un 
volume de donnees sismiques traitees, oil la quantite associee a chaque voxel dans le 
volume de donnes traitees est le resultat obtenu par Texecution de l'operation de 
lissage dans le voxel dans le volume de donnees sismiques. 

Le document WO 02/13139 Al decrit un procede de traitement damages 
20 sismiques comprenant les etapes : obtenir un ensemble de donnees d f une image 
initiale bidimensionnelle ou tridimensionnelle, dans lequel chaque element de 
Tensemble de donnees est 1'intensite initiale d'image du point de 1'image ; calculer 
pour chaque point les derivees partielles de reiement dans n directions pour obtenir 
un ensemble de valeurs derive des derivees partielles; calculer pour chaque point une 
25 matrice carree structurale symetrique k partir des valeurs des derivees partielles; 
effectuer une iteration point h point ponderee par une variable voisine de zero lorsque 
le point est voisin d'un bord et voisine de 1 lorsqu'il est trfes eioigne d'un bord ; et 
repeter ces etapes un nombre de fois pour obtenir l'image traitee. 

Un premier but de Tinvention est d ! ameiiorer le lissage des proprietes de 
30 subsurface dans une structure geologique representee par des mesures sismiques. 



Un deuxteme but de 1'invention est de permettre un lissage simple et rapide 
des propridtes de subsurface, tout en ne lissant pas les discontinuites. 

L!inve^qn_a^ subsurface 
dans une structure gdologique representee par mesures sismiques, dans lequel on 
5 construit une fonction continue S ijfk (t) par interpolation ou approximation des traces 
sismiques discretes d ! une matrice sismique multidimensionnelle, ladite fonction etant 
d6sign<Se comme "trace sismique locale continue", comportant les Stapes suivantes : 

a) , utiliser comme ddcalage optimal hij >pq , k de deux traces sismiques locales 
continues voisines S^t) et S pqik (t), la valeur de d<Scalage rendant maximale leur 

10 fonction de correlation ; 

b) retenir comme voisinage conditionnel d'une trace sismique locale continue 
"centrale" Sy, k (t) le sous-voisinage consistant en des traces voisines Spq, k (t) 
correspondant k des decalages optimaux h% m x associes a des correlations 
Rij.pq^Chij.pq.k) sup^rieures k un seuil predetermine compris entre 0 et 1 ; 

15 c). choisir une propriete de subsurface a lisser au voisinage conditionnel d*un 

point (i, j, k) d'une trace sismique locale continue "centrale" de reference ; 

d) . d£caler les proprietes de subsurface du voisinage conditionnel en 
translatant la variable courante de la valeur du d&alage optimal hy,pq fk ; 

e) . prendre comme valeur liss^e au point (i, j, k) une moyenne des propri&£s 
20 de subsurface decal6es k Tetape d), 

Selon d ! autres caract6ristiques de invention : 

- la moyenne des propri&ds de surface de l'6tape e) est une moyenne qui peut 
Stre pond£ree, par exemple par la valeur de correlation maximale correspondant au 
d€calage optimal, 

25 - la moyenne des propri&ds de subsurface de l'&ape e) est choisie dans 

Tensemble suivant : moyenne arithmetique, moyenne g6om<Strique, moyenne 
harmonique, 

- on choisit comme propriety de subsurface k lisser Tamplitude r£ftechie et 
capttSe par des g£ophones, 
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- si necessaire, Pensemble du processus de lissage est r£p£te un certain 
nombre de fois pour ameliorer le lissage, 

- on visualise des coupes de matrice multidimensionnelle de propriety liss£e 
sur un 6cran de visualisation. 

5 L'invention a egalement pour objet un dispositif pour la mise en oeuvre du 

proced£ selon 1'invention, comportant des moyens pour utiliser corame decalage 
optimal de deux traces sismiques locales, continues voisines la valeur de decalage 
rendant maximale leur fonction de correlation, des moyens pour ddcaler les 
propri€t6s de subsurface du voisinage conditionnel en translatant la variable courante 

10 de la valeur du decalage optimal hy fPq ,k s des moyens pour choisir une propriete de 
subsurface & lisser au voisinage conditionnel d'un point (i, j, k) d'une trace sismique 
locale continue "centrale" de reference, des moyens pour d^caler les proprietds de 
subsurface du voisinage conditionnel en translatant la variable courante de la valeur 
du decalage optimal hij.pq.ket des moyens pour prendre comme valeur lissee au point 

15 (i, j, k) une moyenne des propriety de subsurface decal£es & l'&ape d). 

Selon d'autres caracteristiques de Tinvention : 

-le dispositif comporte des moyens de memorisation et des moyens de 
visualisation de paramfctres sismiques determines & Taide du precede selon 
Tinvention. 

20 Uinvention a enfin pour objet un produit-programme d'ordinateur, 

comportant des elements de code de programme pour ex£cuter les Stapes d'un 
proc6de selon Tinvention, lorsque ledit programme est execute par un ordinateur. 

L'invention sera mieux comprise grSce h la description qui va suivre donn£e k 
titre d'exemple non limitatif en reference aux dessins annexes dans lesquels : 

25 - La figure 1 reprdsente schematiquement une matrice sismique 

tridimensionnelle et une trace sismique locale continue. 

- La figure 2 repr^sente sch€matiquement, un exemple de voisinage local Ny >k 
consistant en un ensemble de trace sismique locale continues S pq ,k(t) voisines de la 
trace centrale de reference Sy, k (t) elle m6me situ£e a la verticale du g£ophone Gy et 

30 interpolant les amplitudes sismiques autour de t=k. 



- La figure 3 represente schematiquement, une section verticale de voisinage 
local conditionnel Nij |k (r) d'une trace sismique locale continue de reference Sij, k (t). 

_ La figure 4 r epresen te scM^tiguement, unej^tip^^ matric e 
sismique tridimensionnelle avec un decalage optimal hy^.k et un voisinage 
5 conditionnel Nij t k(r). 

- La figure 5 represente schematiquement, une section verticale d'une matrice 
sismique tridimensionnelle analogue k la figure 4 avec presence d'une faille ou 
discontinuity 

La figure 6 represente schematiquement, un organigramme fonctionnel d'un 
10 procede selon l'invention. 

En reference aux figures 1 & 5, les elements identiques ou fonctionnellement 
equivalents sont designes ou reperes de maniere identique. 

Sur la figure 1, une matrice sismique tridimensionnelle est obtenue par releve 
de mesures enregistrees par des geophones Gy disposes sur un reseau x,y en des 
15 points de coordonnees i, j. Le releve des mesures echantillonnees dans le temps est 
represente suivant un axe t descendant representatif de la profondeur ou d r une 
verticale descendante k partir de la surface de la terre ou de la mer. Les mesures sont 
caracterisees par leur amplitude, par exemple une amplitude relevee par le geophone 
au temps ou k la profondeur d'echantillonnage t k . La mesure discrete effectuee par 
20 le geophone Gy au temps ou k la profondeur t k est appeiee amplitude sismique Sjj, k . 

Uensemble des amplitudes sismiques correspondant k un geophone Gy de 
coordonnees i, j est une matrice unidimensionnelle (Syi, Sij2,...., Sijk> — , Syw) 
appelees trace sismique discrete, car cette matrice unidimensionnelle correspond k la 
trace selon le point de coordonnees horizontals i, j de la matrice sismique 
25 tridimensionnelle obtenue par mesures sismiques. 

L'axe vertical t oriente selon une verticale descendante designe usuellement le 
temps, mais peut egalement §tre traite pour representer une profondeur k partir de la 
surface. 

L'invention concerne aussi bien 1'application k une troisifcme coordonnee t 
30 representative du temps, qa'k une troisifeme coordonnee t representative de la 
profondeur. 



A partir de la trace sismique discrete situde & la verticale d'un g£ophone Gy 
on ddfinit, par interpolation ou approximation des valeurs discretes autour de t=tk=k, 
une fonction continue Sy^Ct) qui est designee comme "trace sismique locale 
continue". Les m&hodes ^approximation ou d'interpolation de valeurs discretes pour 
5 engendrer une fonction continue sont nombreuses, et comprennent notamment les 
interpolations ou approximations polynomials, ainsi que les interpolations ou 
approximations trigonom&riques polynomiales. 

Toute autre variante d'interpolation ou ^approximation fournissant une 
fonction continue peut £galement Stre appliqu^e a la presente invention pour fournir 
10 une "trace sismique locale continue". 

Sur la figure 2, plusieurs traces sismiques locales continues d^finissent un 
exemple de voisinage d'une trace sismique locale continue "centrale" Sij fk (t) de 
reference. Le voisinage d'une trace sismique locale continue Sy,k(t) est ddfini comme 
Tensemble des traces sismiques locales continues dont les indices spatiaux 
15 horizontaux p,q sont voisins des indices spatiaux horizontaux i, j de la trace sismique 
locale continue de reference. 

A titre d'exemple, les coordonndes spatiales horizontales p,q correspondant au 
geophone G pq sont voisines des coordonndes spatiales horizontales i, j correspondant 
au geophone Gy si les valeurs absolues des differences i-p et j-q sont inferieures & des 
20 entiers donnas, par exemple a 2. 

Dans ce cas, comme repr£sente k la figure 2, la trace sismique locale continue 
Sij,k(t) est associ6e h huit traces sismiques locales continues voisines entourant la 
trace sismique locale continue "centrale" Sy,ic(t). 

Dans le cas de traces sismiques locales continues produites par des mesures 
25 sismiques, le profil des horizons geologiques introduit des d<Scalages verticaux entre 
des traces sismiques locales continues voisines. Dans le but de determiner des 
relations entre deux traces sismiques continues locales voisines Sy,k(t) et S pq>k (t) 
centres sur la m6me coordonn£e verticale d'£chantillonnage t=t k =k et correspondant 
h des coordonn^es spatiales i, j et p,q diff&entes, on calcule la fonction de correlation 
30 Ry^icCh) des deux traces sismiques locales continues. 



La fonction de correlation de deux traces sismiques locales continues 
voisines est obtenue par la formule suivante 



^-..-.-^y.^- - 

oO le numerateur correspond k la fonction de covariance de Si^OO et Spq,k(t) 

tk+A 

obtenue par l'expression suivante Cij,pq^(h) = Js ij k (t) . S wk (t + h) . dt . 



tk-A 



Dans cette integrate definissant Ci jtPq ,k(h) , le parametre A definit une «fenetre 
verticale d' investigation » autour de t=tk==k 

Par exemple, si les traces sismiques locales continues Sij, k (t) et S pq>k (t) sont 
10 des polyndmes trigonometriques de la forme suivante interpolant les donnas 
sismiques 

Sij.k(t) = X a ^ xos ( s ^ t ) +b 3 ,k - sin ( sfi)t ) 

s=l 

S«.k(t) = Xaf .cos(s<Dt) + b** .sin(sa)t) 



s-\ 



15 alors on peut montrer que la fonction de covariance C^ pqi ]c(h) est elle mSme un 
polyn6me trigonometrique de la forme suivante oil les coefficients A s 1J%pqrk et B s li,vqM 
dependent des coefficients a.** a,*** W iM et b,"* : 

Qj,pqjc(h) = XAl j ^\cos(s6)h) + B^ k .sin(s6)h) 

20 Un resultat mathematique classique et connu est que la fonction de correlation 

R^vOi) = " ,pq traduit une similarity des fonctions cornStees 

lorsque cette fonction de correlation se rapproche de L 

L'etude des fonctions de correlation des traces sismiques locales continues 
25 Si jtk (t) et Spq, k (t) permet de definir comme d£calage optimal la valeur hy,^ de h 
correspondant au maximum de la fonction de correlation %pq,k(h) , c'est-k-dire au 
maximum de la fonction de correlation le plus proche de L Contrairement & Fart 
existant base sur une formulation discr&te de la fonction de correlation Rij tPq ,k(h), 




Futilisation d'une formulation continue des traces sismiques et done de la fonction 
de correlation permet d'obtenir un decalage optimum qui n'est pas contraint a Stre un 
multiple entier du pas d'echantillonnage suivant Taxe vertical correspondant a la 
variable t. 

5 La definition du decalage optimal hy tPq ,k de deux traces sismiques locales 

continues voisines permet d ! obtenir une premifere approximation des horizons passant 
par cette trace sismique locale continue de reference Sy,k(t), comme repr^sente aux 
figures 3 k 5. Contrairement a Tart existant, le fait d'utiliser des decalages optimaux 
qui ne sont pas contraint k etre un multiple entier du pas d'echantillonnage suivant 
10 Taxe vertical permet d'eviter des erreurs connues sous le nom « d' aliasing ». 

Sur ces figures 3 k 5, on definit un voisinage conditionnel Ny^Cr) de la trace 
sismique continue locale de reference Sy,k(t) qui est un sous-voisinage du voisinage 
initial des traces sismiques locales continues retenues pour le calcul de correlation et 
de definition des decalages optimaux. Le voisinage conditionnel NyjcOO est choisi de 
15 maniere que, pour toute trace sismique continue locale S pq> k(0 appartenant k Ny,k(r), 
le decalage optimal hy fP qjc correspondant est tel que la correlation Rij,pq,k(hy,pq,k) entre 
Sy ? k(t) et Spq,k(t) est superieure k un seuil predetermine r compris entre 0 et 1. 

Comme on le voit sur la figure 5, la fixation du seuil r permet egalement 
d'englober dans un voisinage conditionnel une discontinue ou une faille geologique, 
20 ce qui constitue un avantage important par rapport k Tart anterieur. 

Les dispositions precitees fournissent ainsi une modeiisation continue 
permettant la mise en oeuvre de l'invention. 

En reference k la figure 6, un procede de lissage ou propriete de subsurface 
dans une structure geologique represente par des mesures sismiques comporte une 
25 premiere etape 100 de modeiisation numerique, comme decrit en reference aux 
figures 1 a 5, pour definir des traces sismiques locales continues d'une matrice 
sismique multidimensionnelle, calculer les decalages optimaux de traces sismiques 
locales continues voisines et definir les voisinages conditionnels des traces sismiques 
locales continues "centrales" de reference. 

30 Apr£s cette premiere etape 100 de modeiisation numerique, une etape 101 

permet de choisir la propriete k lisser sur un voisinage conditionnel. Cette propriete k 




lisser sur un voisinage conditionnel peut consister en toute fonction 
multidimensionnelle de la subsurface k analyser. 

En particulier, la propri6te k lisser peut simplement consister^n 1' amplitude 
du signal sismique reflechi par les horizons de la subsurface en direction des 
g6ophones de mesure. 

A l'etape 102, on remplace les propriety du voisinage conditionnel par les 
propridtes decaltSes de la valeur du decalage optimal h^k, de manifcre k constituer 
un ensemble de propriety relatives au meme horizon. 

Ainsi, dans le cas d"une trace locale sismique continue centrale de reference 
Sij, k (t), on remplace toutes les valeurs Spq tk (t) du voisinage conditionnel par les 
valeurs Spq,k(t+hij,pq,k) - 

On effectue ensuite k l'etape 103 une moyenne des valeurs decalees 
Spq^t+hij,^^ etablies k l'etape 102 et on choisit comme valeur de la propri&e au 
point de coordonnees spatiales i,j et de coordonnee temporelle ou de profondeur 
t=tic=k , la moyenne ainsi calcuiee. 

^invention s'applique k plusieurs variantes d'etablissement de moyennes : on 
peut par exemple utiliser une moyenne pond£r£e ou non ponderee; ou, 
alternativement, on peut egalement utiliser une moyenne arithmetique, une moyenne 
g£om€trique, une moyenne harmonique , ponder£e ou non pond6r6e. 

Dans le cas d'une moyenne ponder<Se, ont peut avantageusement utiliser 
comme f acteur de pond€ration la valeur maximale de la fonction de correlation 

Rij,pq,k (hij,pqjc) , en refietant ainsi le degre de confiance accordable k chaque 
decalage optimal conrespondant k un maximum d'une fonction de correlation. 

Le precede de lissage ou de filtrage d'une propriety selon la prdsente 
invention pennet ainsi d'ameliorer le contraste et la precision des proprietes de 
subsurface dans une matrice sismique multidimensionnelle. 

Si necessaire, afm d'ameliorer le lissage, l'ensemble du processus ddcrit par 
la presente invention peut Stre repete sur le lissage obtenu lors d'une iteration 
preeddente. L'etape 104 cOnstitue une etape iterative dans laquelle l'operateur peut 
definir le nombre derations adequates k la resolution de son probteme. 




L'invention est de preference mise en ceuvre sur un dispositif comportant des 
moyens pour accomplir successivement les etapes du proc€d6 ddcrit en reference k la 
figure 6. 

En particulier, ce dispositif comporte des moyens de visualisation des 
5 parametres sismiques ou proprietes obtenues & Tetape 105 et des moyens de 
memorisation de programme et de calculs intermediaires* 

De preference, ce dispositif est un dispositif programmable commande par un 
produit-programme d'ordinateur, realise pour mettre en ceuvre un procede selon 
l'invention. 

10 L'invention decrite en reference h plusieurs objets particuliers n'y est 

nullement limitee mais couvre au contraire toute modification de forme et toute 
variante de realisation dans le cadre et l'esprit de invention, Tessentiel etant d'utiliser 
des valeurs de proprietes recalls par rapport aux decalages optimaux determines 
comme valeurs correspondant S. la maximisation des fonctions de correlation de 

15 traces sismiques locales continues ou mSme discr&es. 




RT?VF.NT>TCATIONS 

1. Procede de lissage d'une propridte de subsurface dans une structure 
ggologique representee par mesures sismiques, 

continue S ijlk (t) par interpolation ou approximation des traces sismiques discretes 
d'une matrice sismique multidimensionnelle, ladite fonction S ijtk (t) etant designee 
comme "trace sismique locale continue", comportant les etapes suivantes : 

a) , utiliser comme decalage optimal de deux traces sismiques locales 
continues voisines S ij>k (t) et Spq^t), la valeur de decalage rendant maximale leur 
fonction de correlation, ce decalage optimal n'etant pas contraint a etre un multiple 
entier du pas d'echantillonnage vertical ; 

b) . retenir comme voisinage conditionnel d'une trace sismique locale 
continue "centrale" Si jjc (t) le sous-voisinage consistant en des traces voisines S pq ^(t) 
correspondant a des decalages optimaux associes a des correlations Ry,pq,k(h) 
superieures a un seuil predetermine compris entre 0 et 1 ; 

c) . choisir une propriete de subsurface a lisser au voisinage conditionnel d'un 
point (i, j, k) d'une trace sismique locale continue "centrale" de reference ; 

d) . decaler les proprietes de subsurface du voisinage conditionnel en 
translatant la variable courante de la valeur du decalage optimal (hjj,pq,k) ; 

e) . prendre comme valeur lissee au point (i, j, k) une moyenne des proprietes 
de subsurface decaiees & l'etape d). 

2. Procede selon la revendication 1, dans lequel la moyenne des proprietes de 
surface de l'etape e) est une moyenne ponderee, par exemple par la valeur de 
correlation maximale correspondant au decalage optimal. 

3. Procede selon la revendication 1 ou 2, dans lequel la moyenne des 
proprietes de subsurface de l'etape e) est choisie dans l'ensemble suivant : moyenne 
arithmetique, moyenne geometrique, moyenne harmonique ponderee ou non. 

4. Procede selon la revendication 1, dans lequel on choisit comme propriete 
de subsurface a lisser, l'amplitude refiechie et captee par des geophones. 

5. Procede selon l'une quelconque des revendications precedentes, dans lequel 
on itere le processus de lissage en l'appliquant sur le resultat d'un lissage anterieur. 



6. Precede selon la revendication 1 ou la revendication 5, dans lequel on 
visualise une matrice multidimensionnelle de propriety liss£e sur un dcran de 
visualisation. 

7. Dispositif pour la mise en oeuvre du procede selon Tune quelconque des 
5 revendications 1 & 4, comportant des moyens pour utiliser comme decalage optimal 

de deux traces sismiques locales continues voisines la valeur hy,pq,k de decalage 
rendant inaximale leur fonction de correlation, des moyens pour decaler les 
propri6tes de subsurface du voisinage conditionnel en translatant la variable courante 
de la valeur du d€calage optimal Ifypq^ , des moyens pour choisir une propria de 
10 subsurface a lisser au voisinage conditionnel d'un point (i, j, k) d'une trace sismique 
locale continue "centrale" de reference, des moyens pour decaler les propriety de 
subsurface du voisinage conditionnel en translatant la variable courante de la valeur 
du d£calage optimal hty.pq^ et des moyens pour prendre comme valeur lissde au point 
(i, j, k) une moyenne des propri6t£s de subsurface ddcalees a I'etape d). 

15 8. Dispositif selon la revendication 7, comportant des moyens de 

memorisation et des moyens de visualisation de paramfetres sismiques determines k 
l'aide du precede selon Tune quelconque des revendications 1 & 6. 

9. Produit-programme d'ordinateur, comportant des elements de code de 
programme pour executer les etapes du procede selon Tune quelconque des 

20 revendications 1 a 5, lorsque ledit programme est execute par un ordinateur. 

10. Produit-programme d'ordinateur, comportant des elements de code de 
programme pour executer les etapes du procede selon la revendication 6, lorsque 
ledit programme est execute par un ordinateur. 
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